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Nis tony Detailed Current Edit areerr 

CLARATIONS ; Declarative Par of odule 

H DAS = Standard Double Precision Floating DASIN 
MTHSDASIN_R7 = Special DASIN routine 
MTHSDASIND = Standard Double Precision Floating DASIND 
MTHSDASIND_R7 = Special DASIND routine 
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; Double Precision Floating Point Sine r 16-SEP-1984 01:14:07 \VAX/VMS Macro v04-00 Page 1 
° aes ter es 14 ERTHRTL. SRETATHDASIN.MARS 1 ’ (1) 
-TITLE MTHSDASIN ; Double Precision Floating Point Sine routine 
3; (DASIN,DASIND) 
IDENT /1-006/ : File: MTHDASIN.MAR Edit: RNH1007 


RERARRAEAEREAAE EAA AAAA AA AAAAEAEAAAAEAARAAAAAAAEAAAEEEEEEEEAEREREAAAAAAAAAAEREREE 


Se COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 


DOoODOOOOCSooO 


OOOCCOOOCoo 


DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
RIGHTS RESERVED. 


a 
4 
5 
7* ® 
| § 3? ® 
3* ® 
33g Me 
PY ® 
000 11 ;* THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
8 1 3* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
13 ;* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 

4 14 ;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
000 15 ;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 

0900 16 ‘* TRANSFERRED. * 

3 ® 

000 18 ;* THE INFORMATION IN THIS SOFTWARE iS SUBJECT TO CHANGE WITHOUT NOTICE ® 
0000 19 ;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT ® 
44 Y :* CORPORATION. * 
3 ® 

0000 § 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
4 ? :* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. * 
3 ® 

0900 i s conan GuadulehisenanineaneenniehianehnepNneadeabakedenedeetabedenénneanenes 

0000 8 ; 

0000 9 ; FACILITY: MATH LIBRARY 

0000 0 ;++ 

464 1 ; ABSTRACT: 

44 : : MTHSDASIN is a function which returns the double precision floating point 
000 4 ; arcsine in radians of its double precision floating point argument. The call 

8308 2 3 iS Standard call-by-reference. 

0000 7 : MTHSDASIND is a function which returns the double precision floating point 
000 3; arcsine in radians of its double precision floating point argument. The call 
000 9 ; is standard call-by-reference. 

000 40 ; 

000 41 ;-- 

000 $s : 

000 re 3; VERSION: 01 

$00 45 ; HISTORY: 

O60 rt 3 AUTHOR: 

00 ? 3 Peter Yuo, 29-Jun-77: Version 01 

09 2 ; MODIFIED BY: 
0 0; 

00 13 

00 e3 


MTHSDASIN : pte Precision Floating P int” sine tr 16-SEP-1984 01:14:07 VAX/VMS Macro v04-00 Page 
1-006 furs? ; Detailed Current Edit History g- set a MTHRTL. SR aT MTHDASIN.MAR; 1 . 3) 


-SBTTL HISTORY ; Detailed Current Edit History 


P 


ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: none 
Edit History for Version 01 of MTHSDASIN 


0-2 +t ie oe A to B  saceetntatas 
eee Chan 
Changed errorshandt ing. mechanism. Put error result in RO:R1 before 
galt ng MTHSS$SIGNAL in order to allow user modif das: r result. 
ix a. involving wrong addressin sete. JMT 9- 
Upde e version number and eoere A ae neree 188) “16-Nov-78 
qgnae.. gis INVARG to MTHSK_INVARGMAT o-iee Re 78 
o the PSECT directive. 
Declare Bin B+ SBL 1/7-May-1 
roange JSB entry to MTHSDASIN R7. RBG 28- ~ pasaalthda 
Added degree entry points. RNH 29-MAR-1981 


COOCCCCOCOCOOOOOOOCoOoOoO 
o 
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ao 
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NS ; Declarat ve Port of Modul aes 7 bere 14 UMTHRTL . SREIATHDASIN.MAR: 1 ’ ) 


-SBTTL DECLARATIONS ; Declarative Part of Module 


=a 
ov 


I 


; program section for math routines 


INCLUDE FILES: MTHJACKET.MAR 
; EXTERNAL SYMBOLS: 
0 -DSABL GBL 
0 -EXTRN MTHSDSQRT_RS 
8 ~EXTRN MTHSDATAN_R7 
-EXTRN MTHSDATAND R7 
00 sEXTRN MTHSK_INVARGMAT 
BOS -EXTRN MTHSSSIGNAL 
00 ; 
88 ; EQUATED SYMBOLS: 
00004080 000 38.1.0 2 *#1.6 : 1.0 
00000004 B08 value =4 : value.rd.r 
000 3 
4 ; MACROS: none 
990 > PSECT DECLARATIONS: 
0000000 .PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
00 


: OWN STORAGE: none 
: CONSTANTS: 
DL 


PI_OVER 
68C2 A221 OFDA 40C9 . 


2: 
GORD *0040311, *0007732, “O1Z108 | 5, “OOSAI0e 
68C2 A221 OFDA 4149 


00000000 00004384 
00000000 00004434 


-LONG  *X00043B4, “x0 : 90 
- "LONG = *X0004434, “x0 : 180 


SCBHVOOCOWMDHMDOOOOCOCOSOOSOSOSOOSSCSSSSOSSOOOOOCOOOOOOOOOOOOO 

3 tt ~QDOODODOOOOOWOO 000000000 090909 0909 09090909 II NINN 

OWONAUE WN OOD NAME WW SOO NA NEW OO ODONAULS WN OOONOULSW 
7 


ed ae a te a tw 


SOoooooooooooooooooo 
SOoooooooooooooooo 


Nee | “OOOO OOOOSoSoSo 


D_PI: 
= «WORD %*0040511, “0007732, “whew *0064302 
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-SBTTL MTHSDASIN = 


eB 


P=-1984 01:14:07 VAX/VMS Macro V04-00 Pa 4 
poms ri a tat tt 


tandard Double Precision Floating DASIN 


oc 


Arhs 


++ 

FUNCTIONAL DESCRIPTION: 

; DASIN = double precision floating point function 
; DASIN(X) is computed as: 


If X = 0, then DASIN(X) = 0. 

If X = 1, then DASIN(X) = PI/2. 

If X = =f, then DASIN(X) = =PI/2. 

If 0 < ix! <1, then DASIN(X) = ATAN(X/SQRT(1-X##2)). 
If 1 < ixXi, error. 


OOOCCCOCOOCOOOOoOoOoOoO 


PMOPOITOPIPOPINIPIPIPINIPOPONIPUNoNonoreornrnornn 


OOOCCCOCCOCCOOOOCOCOCOOOOOOOO 


; CALLING SEQUENCE: 
DASIN.wd.v = MTHSDASIN(x.rd.r) 
INPUT PARAMETERS: 


LONG = 4 3; define senguare multiplier 
x = 1 * LONG ; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: double precision floating arcsine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: MTHS_INVARG if {Xi > 1 with reserved operand in RO/R1 
(copied to the signal mechanism vector CHFSL_MCH_RO/R1 by LIBSSIGNAL). 


Associated message is: ‘INVALID ARGUMENT’. Resul€ is reserved 
operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 


00000004 
00000004 


SSssssssss SDOOQOCOOCOSOSOOOOOOOOOOOSooOOO 
UNO ODNA UNE WO OD NAULS WN OOO NAULS WN OUOONOUS WN 7H 


POPOIPOPIPIPIPININININPININPINININININININY 


DWODOMMMNNNOOCOCOCCOOCCOOCOOCOOCOOCSOCOCoOooOooo 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


40FC ENTRY MTHSDASIN, “M<IV, R2, R3, R4, R5, Rb, R7> 
3 standard call-by-reference entry 
3; disable DV (and FU), enable IV 


MTHSFLAG_JACKET ; flag that this is a jacket procedure in 
MOVAB G*MTHSSJACKET_HND, (FP) 


ce ee ee ee ed ed ed od dd dd wd dd 3 ed 3 dd 2d 3 dd 3 


m—O DONO 


6D 00000000°GF 9E 


set handler address to jacket 
handler 


SOOooooooooooooooooa 


SoOOoOOooOooOooOoCoooooooSo 
POPSPOROPORNONOPoPUnonononys 


vs 


miggesi a SOME THOE Hts F ERERIBEL CUDHSSG WARE MESRBIGBR aes PO 


: igre of an error in reyes a 
: iF feats convert signal to user PC 
50 04 BC 70 176 MOVD — avalue(AP), RO : ROURI = {Rt = avalue(AP) 
| 17 BSBB  =s- MTHSDASIN_R7 ; call special DASIN routine 
04 + f V8 RET 3; return with result in RO/R1 


a2. Ses. Weanie 
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1-006 ATH DASIN_R? = Special DASIN routine ~SEP=1984 1:34:93 MTHRTL.SRCIMTHDASIN.MAR; 1 . (8) 
} 1 ~SBTTL MTHSDASIN_R7 = Special DASIN routine 
! $ ; Special DASIN = used by the standard routine and direct JSB call. 
185 ; CALLING SEQUENCE: 
1 § : save anything needed in RO:R7 
3 1 3 MOVD eo ; input in RO/R1 
1 § ; JSB MTHSDASIN_R7 
8 ! 3 RSB 3; return with result in RO/R1 
0 191 ° 
136 MTHSDASIN_R7:: 3 special DASIN routine 
195 MTHSDASIN_R9I:: ; Release 1 name 
56 650 9 0 9 194 MOVD RO, R6 3 save X in R6/R7 
2 13 00 195 BEQL RETURN : return DASIN(O) = 0 
3s 5 139 
0 197 ; 
00 198 ; O< Ix! 
8 199 ; 
035 00 
50 8000 8F AA tH 01 BICW #*x8000, RO 3; RO/R1 = ix! 
08 50 71 00 A 0¢ CMPD RO, S*#SD_1.0 : compare !X! with 1.0 
17 18 0030 20 BGEQ GEd_T0_1.0 : branch if {X! >= 1.0 
Ser 04 
003F 05 ; 
003F 06; 0< Ix! < 1.0 
003F 07 ; 
003F 08 
50 50 64 0O03F 09 MULD RO, RO 3 RO/R1 = X#e2 
50 08 50 63 0042 10 SUBD RO, S$*#SD_1.0, RO 3 RO/R1 = 1.0 = X#e2 
00000000 "EF 16 0046 11 JSB MTHSDSQRT RS 3; RO/R1 = DSQRT(1-X**2) 
0 6 50 67 004C \¢ DIVD3 RO, R R 3 RO/R1 = X/DSQRT(1-X*#2) 
00000000‘ EF 17 boas \7 JMP MTHSDATAN_R7 3: RO/R1 = DATAN(X/DSQRT(1-X**2)) 
0056 15 
Baee 16 ; 
056 173 1 s¢ ix! 
0056 18 3 
B20 1 
56 0 GEQ_T0_1.0: 
oc 614 B28 1 BGTR ERROR : branch to ERROR if {Xi > 1.0 
05 j 
O58 958 | 
° 5 : ix! = 1.0 
33 559 
50 AS AF A 5 8 MOVD D_PI_OVER_2, RO 3: RO/R1 = PI/2 
36 73 «005C 9 TSTD RB ; test the sign of x 
2 1 5 0 GEQ my une 3: branch if X > 
50 S fF 6 1 MNEGD RO, RO 3: RO/R1 = -PI/2 
0 o7 § RETURN: RSB 3; return with result in RO/R1 
4; 
4 5; 1< txt, error 
ae 
4 


M14 
MTHSDASIN poe ble Precision Floating Point Sine r 1 “$ oer 138% 91: }: 97 AX/VMS Macro Vv04-00 Page 7 
1 - DASIN_R? = Special DASIN routine 4 9 (CMTHRTL.SRCJMTHDASIN.MAR; 1 (5) 


E 4 ERROR: PUSH (SP) : retuen. PC from JSB routine 
E 00° SF 6 Oey M5 iL erenee INVARGMAT, -(SP) : io" tion value 
01 F A #1 T, RO : = = cooule =. = reserved operand -0.0 

E 4 5 goes gnal_ mechanism vector 
E 4 ; Cnese MCH_R PRI) so error handler 
E 4 : can modify the result. 

00000000'GF 02 44 CALLS #2, G*MTHSSSIGNAL : $ “—e error and use real user's PC 

45 ; independent of CALL vs J 
ee $$ RSB ; return - RO restored from CHFSL_MCH_RO/R1 
076 48 
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"nn 
—@ 
ov 


=006 hTHSDASIN 


? -SBTTL MTHSDASIND = Standard Double Precision Floating DASIND 
, p++ 
i ¢ : FUNCTIONAL DESCRIPTION: 
1 2$ 3 DASIND = double precision floating point function 
O76 28 ; DASIND(X) is computed as: 
076 260 ; If X = 0, then DASIND(X) = 0, 
13 ] 3 If X= then DASIND(X) = 90, 
? ¢ 3 If X = ={, then DASIND(X) = =90. 
076 65 ; If 0 < iXi < 1, then DASIND(X) = ATAND(X/SQRT(1-X*##2)). 
Baee ¢: : If 1 < iXi, error. 
0076 66 ; CALLING SEQUENCE: 
Boe re: ; dasind.wd.v = MTHSDASIND(x.rd.r) 
0076 270 : INPUT PARAMETERS: 
00000004 0076 % LONG = 4 ; define Lonquere multiplier 
00000004 onee te x = 1 * LONG ; Contents of x is the argument 
Boe a : IMPLICIT INPUTS: none 
0076 377 : OUTPUT PARAMETERS: 
0076 4s : VALUE: double precision floating arcsine of the argument 
One ei ; IMPLICIT OUTPUTS: none 
0% es ; COMPLETION CODES: none 
O76 285 ; SIDE EFFECTS: 
0076 89 : Signals: MTHS_INVARG if {X! > 1 with reserved sporend in RO/R1 
0076 88 ; (col80ed to the signal mechanism vector CHFS$L_MCH_RO/R1 by LIBSSIGNAL). 
076 89 ; Associated message is: ‘INVALID ARGUMENT’. Result is reserved 
gfe 39 ; operand -0.0 unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1 
0076 5 ; NOTE: This procedure disables floating point underflow, enables integer 
143 93 ; overflow. 
7 94 ; 
076 99 32° 
b76 599 
40FC 007 98 -ENTRY MTHSDASIND, “M<IV, R2, R3, R4, RS, Rb, R7> 
? 99 ; standard call-by-reference entry 
7 8 : disable DV (and FU), enabl 
e 1 MTHSFLAG_JACKET ; flag that this is a jacket procedure in 
6D 00000000'GF 9€ 7? MOVAB G*MTHSSJACKET_HND, (FP) 
grr ; set handler address to jacket 
4 ; handler 


ee 


riggs ARRAN NSLLaM SOP ALS CRETE NIE PRAM. SRESRABHEAR nas PG 


case of an error in retin 


7F Py 
7F § ; If an error convert signal to user PC 
4 3; and resigna 
50 04 6C 70 7F 5 MOVD @value(AP), RO 3 RO/R1 = {Xi = @value(AP) 
01 10 0083 306 BSBB MTHSDASIND. R7 3 call special DASIND routine 
04 0085 307 RET 3; Return with result in RO/R1 


— = tee 
po a | el 
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MTHSDASIN ; Do r 
1-006 = Special DASIND routine “-SEP-19 MTHRTL.SRCJMTHDASIN.MAR; 1 ( 


Pr 
D 


° 19 -SSTTL MTHSDASIND_R7 = Special DASIND routine 
6 1 ; Special DASIND = used by the standard routine and direct JSB call. 
14 : CALLING SEQUENCE: 
6 2 3 save anything needed in RO:R7 
086 16 j MOVD Py 3; input in RO/R1 
086 317; JSB MTHSDASIND_R7 
+4 18 ; RSB ; Return with result in RO/R1 
ae 
$38 MTHSDAS!ND_R7:: 3 special DASIND routine 
56 650 19 0 MOTD RO, R6 3 save X in R6/R7 
2F 1 089 BEQL D_RETURN ; Return DASIND(O) = 0 
0088 355 ; 
0088 6:0 < Ix! 
BB 
50 8000 8F AA 0088 23 BICW #*x8000, RO ; RO/R1 = IX! 
08 50 71 0090 0 CMPD RO, S*#SD_1.0 3 compare iX: with 1.0 
17 18 0093 «= (331 BGEQ 0_GEQ_T0_T.0 > branch if !X!} >= 1.0 
$33 
0095 334: 0 < Ix! < 1.0 
$3: 
50 50 64 0095 37 MULD2 RO, RO ; RO/R1 = X**#2 
50 08 50 63 0098 335 SUBDS RO, S*#SD_1.0, RO 3 RO/R1 = 1.0 = Xe*2 
QOOOOOOO'EF 16 009C¢ 39 JSB MTHSDSQRT_RS 3 RO/R1 = DSQRT(1-X**2) 
0 56 SO 67 OQOA2 379 DIVD3 RO, R6, R 3 RO/R1 = X/DSQRT(1-X*#2) 
00000000" EF 17 erat 30 JMP MTHSDATAND_R7 3 RO/R1 = DATAND(X/DSQRT(1-X**#2)) 
OOAC 305 
QOAC 344 ; 
OAC 345 ; 1 =< {XxX} 
Ooar 849° 
QOAC 348 D_GEQ_TO 1.0: 
B6 14 OQOAC 49 BGTR ERROR 3 branch to ERROR if !X! > 1.0 
OAE 50 
OAE 51 
OAE 26 5 
AE 55 3: ix! = 1.0 
at 
50 FFSE gf 7 A 2 MOVD 0.90, RO ; RO/R1 = 90 
Gh i GEE Bein ride 
: branc 
50 5 7 oer 33 , tue RO, RO ; RO/R1 = -90 
05 OBA 61 7 RSB ; Return with result in RO/R1 
53 ie 
pee 64 END 


OE 
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Symbol table . SrSEp-19R6 $1:3ti99 PATUNTS eacsartipAssyemar:1 (2% {}, 
D_180 0000018 R 1 
D-90 18 & 1 
D_GEQ_T0_1.0 AC R 1 
PI 05 R 1 
D_P1_OVER_2 8 1 
D-RETURN OBA R 01 
OR 9 R 01 
GEQ_10_1.0 0000056 R 01 
= 00000004 
MTHSSJACKET_HND teeeeree x (1 
MINSDAgIN © 00000020 RG 6 
NTHISDAS IND Goug00% RG a 
MTHSDASIND_R7 00000 8 RG 1 
MTHSDASIN_R7 00000030 RG 01 
MTHSDASIN RO 00000030 RG 3} 
MTHSDATAND_R7 etereeee YX 0 
MTHSDATAN_R7 eeereeee =X 00 
MTHSDSQRT_RS rerereee x 00 
MTHSK_INVARGMAT eeeereer XX 00 
RETURN 00000063 R 01 
0.1. = 00004080 
vACUE = 00000004 
em ge ee eR ee ee ee ee + 
: Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
WS 00000000 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 000000BB ( 187.) O1¢ 1.) PIC USR CON’ REL LCL SHR EXE RD NOWRT NOVEC LONG 
ee 
! Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 34 68 00:00.09 68 00:01.08 
Command processing 116 Fhe Sh B00 289 «88 
Pass 1 96 00:00:00.99 00:00:03.39 
Symbol table sort 9 0:00:00.01 0:00:00.01 
Pass 2 7 :00:00. ; 8: : 03 
Symbol iable output 3 :00:00. :00:00. 
pases synopsis ge : 9: 3. j ‘ 88 bs 
ross-reference outpu : ® 
Assembler run totals 32 80; 0:02.6 0 80 93-92 


The working set Limit was 750 pages. 
5305 bytes (11 pages) of virtual memory were used to buffer the intermediate code. 


There were 10 pages of symbol table space allocated to hold 24 non-local and 0 local symbols. 
424 source lines were read in Pass 1, producing 14 object records in Pass 2. 
1 page of virtual memory was used to define 1 macro. 
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Macro Library name Macros defined 
_S255$DUA28: (SYSLIBISTARLET .MLB;2 nena a 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDASIN/OBJ=0BJ$:MTHDASIN MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 


AH-BT13A-SE EQUIPMENT CORPORATION 
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